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ABSTRACT 



Context. The metallicity dependence of the Cepheid Pi-relation is of importance in establishing the extra-galactic distance scale. 
Aims. The aim of this paper is to investigate the metallicity dependence of the PL-relation in V and K based on a sample of 68 
Galactic Cepheids with individual Baade-Wesselink distances (some of the stars also have an HST-based parallax) and individually 
determined metallicities from high-resolution spectroscopy. 

Methods. Literature values of the V^-band, A'-band and radial velocity data have been collected for a sample of 68 classical cepheids 
that have their metallicity determined in the literature from high-resolution spectroscopy. Based on a (V^ — K) surface-brightness 
relation and a projection factor derived in a previous paper, distances have been derived from a Baade-Wesselink analysis. PL- and 
PLZ-relations in V and K are derived. 

Results. The effect of the adopted dependence of the projection factor on period is investigated. The change from a constant p-factor 
to one recently suggested in the literature with a mild dependence on log P results in a less steep slope by 0. 1 unit, which is about the 
1-sigma error bar in the slope itself. The observed slope in the PL-relation in V in the LMC agrees with both hypotheses. In K the 
difference between the Galactic and LMC slope is larger and would favour a mild period dependence of the p-factor. The dependence 
on metallicity in V and K is found to be marginal, and independent of the choice of p-factor on period. This result is severely limited 
by the small range in metallicity covered by the Galactic Cepheids. 

Key words. Stars: distances - Cepheids - distance scale 



1. Introduction 

Obtaining accurate distances to stars is a non-trivial matter 
Cepheids are considered an important standard candle as they 
are bright and are thus the link between the distance scale in 
the nearby universe and that further out via those galaxies that 
contain both Cepheids and SNIa. 

Distances to local cepheids may be obtained via main- 
sequence fitting for those Cepheids in clusters (e.g. Feast 1999) 
or via determination of the parallax. Until recently only Polaris 
had an accurate parallax determination via Hipparcos. Benedict 
et al. (2007) recently published absolute trigonometric paral- 
laxes for nine Galactic Cepheids using the Fine Guidance Sensor 
on board the Hubble Space Telescope, and revised, more ac- 
curate, Hipparcos parallaxes have also become available (van 
Leeuwen et al. 2007, van Leeuwen 2007). 

In addition, distances to Cepheids can be obtained from the 
Baade-Wesselink method. This method relies on the availability 
of surface-brightness (SB) relations to link variations in colour 
to variations in angular diameters and understanding of the pro- 
jection (p-) factor that links radial velocity to pulsational velocity 
variations. 

The SB relations can be obtained using Cepheids that have 
accurate interferometrically determined angular diameters over 
the pulsation phase and accurate multi-wavelength light curves. 
This has allowed Kervella et al. (2004) and Groenewegen (2007; 
hereafter G07) to derive such relations, in particular using the 
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{V — K) colour, which gives the highest precision in the derived 
angular diameter 

The p-factor can also be obtained from Cepheids with inter- 
ferometrically determined angular diameters over the pulsation 
phase when in addition accurate RV curves are available (see e.g. 
Merand et al. 2005). G07 assembled all currently available data 
to find that a constant p-factor represents the best fit to the avail- 
able data. A strong dependence on period, like p ^ —0.15 log P 
as proposed by Gieren et al. (2005), could be ruled out. A mod- 
erate dependence of p ~ —0.03 logP as used in the Gieren et al. 
papers (1993, 1997, 1998, Storm et al. 2004, Barnes et al. 2003), 
orp ^ —0.064 log P as advocated by Nardetto et al. (2007) are 
also consistent with the available data. 

Another aspect of the Cepheid PL-relation that is still a mat- 
ter of debate is the metallicity dependence of the Cepheid PL- 
relation, and its dependence on wavelength. Observations seem 
to consistently indicate that metal-rich Cepheids are brighter, 
and various estimates have been given in the literature, —0.88 ± 
0.16 mag/dex (BRI bands, Gould 1994), -Q.A4+^^\ mag/dex 
(VR bands, Sasselov et al. 1997), -0.24 ± 0.16 mag/dex {VI 
bands, Kochanek 1997), -0.14 ± 0.14 mag/dex {VI bands, 
Kennicutt et al. 1998), -0.21 ± 0.19 in V, -0.29 ± 0.19 in 
W, -0.23 ± 0.19 in /, -0.21 ± 0.19 mag/dex in K (Storm et 
al. 2004), and -0.29 ± 0.10 mag/dex {BVI bands, Macri et al. 
2006). 

The potential caveat is that no individual abundance deter- 
minations of individual Cepheids are being used in these studies 
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but rather abundances of nearby Hll regions, or even a mean 
abundance of the entire galaxy. 

Groenewegen et al. (2004) tried to attack the problem by us- 
ing Cepheids with individually determined metallicities, 37 in 
our Galaxy, 10 in the LMC and 6 in the SMC. Considering the 
Galactic sample, a metallicity effect was found in the zero point 
in the VIWK PL-relation of -0.6±0.4 or -0.8±0.3 mag/dex, 
depending on the in- or exclusion of one object. The distances to 
the galactic Cepheids mostly came from BW-distances available 
at that time from Fouque et al. (2003). 

In the last few years significant progress has been made in 
several areas. On the one hand more Cepheids have interferomet- 
rically determined angular diameters. Combined with improved 
theoretical studies (Nardetto et al. 2007, and references therein) 
this has led to an improved understanding of the p-factor and SB 
relations, as was studied in G07. On the other hand, more basic 
photometric light curves and radial velocity data have become 
available and new studies have presented metallicity determina- 
tions for significant samples (see references below). 

In this paper, Baade-Wesselink distances to Cepheids with 
metallicity determinations are revisited, using the p-factor and 
the SB {V — K) colour relation derived in G07, updated as 
described in Sect. 3. Section 2 describes the selection of the 
photometric and radial velocity data. Section 3 outlines how the 
data are modelled. Section 4 describes how the binary Cepheids 
are treated, and new and updated orbital elements are presented. 
Section 5 describes the results. Section 6 presents conclusions. 



2. The sample 

In a first step Cepheids with recent (post-2000) accurate metal- 
licity determinations were selected (from Andrievsky et al. 
2002abc, 2003, 2004, 2005, Luck et al. 2003, 2006, Luck 
& Andrievsky 2004, Kovtyukh et al. 2005a,b, Yong et al. 
2006, Lemasle et al. 2007, Mottini et al. 2008), and for those, 
the McMaster Cepheid Photometry and Radial Velocity Data 
Archivqj was checked for the availability of V, K and radial 
velocity (RV) data. For the stars selected, the literature was 
searched for additional data. 

Table [T] lists the 68 stars selected with the source of data 
indicated. The last column indicates if particular datasets and/or 
ranges in Julian date (JD) were not used in the fitting described in 
the next section to derive the distance. This is typically because 
these datasets show very large residuals, or the stars show a no- 
ticeable change or shift in period, and the range in JD had to be 
restricted, typically centred on the available A'-magnitude data. 
For the determination of the orbital parameters of the Cepheids 
in binaries (Sect|4|i, all available RV data were used. 

For four Cepheids new high-resolution spectra were obtained 
in order to derive the radial velocitjQ Observations took place in 
July and August 2007 using the Coralie spectrograph at the 1 .2m 
Euler telescope located at the La Silla observatory, Chile. This is 
a fibre-fed echelle spectrograph (2" fibres with one on the object 
and one on the sky), which covers the 3880 to 68IOA region in 
68 orders with a spectral resolution of about 50 000. Exposure 
times are between 3 and 10 minutes. During each night, several 
exposures with a tungsten lamp were taken to measure the rel- 
ative pixel sensitivity variation of the CCD. A Thorium-Argon 
lamp was observed at the beginning and end of the night for 



wavelength calibration. The Heliocentric velocities are obtained 
as part of the standard pipeline reduction whereby the spectrum 
is cross-correlated against a mask, in this case that of a G2 star. 
The formal precision on the RVs is of the order 5-15 m/s. The 
RVs are listed in Table|2l 

Like in G07, in a first step to attempt to homogenise the 
datasets, the light curves are read in the program Period04 (Lenz 
& Breger 2005), which allows for an easy visualisation of differ- 
ent datasets. Period04 is used to calculate the period, and phase 
the data to see if there are clear outliers in the data (which are 
flagged and excluded) or clear off-sets between datasets. Some 
datasets were in the end not considered because of very large 
error bars. The rms in the fit was determined for each dataset 
separately and this was used as the typical error bar (unless the 
rms was consistent with the original error bars quoted, in which 
case these were retained). 

X-band data on the CIT/CTIO system were transformed to 
the SAAO-Carter system using the formula in Carter (1990), 
and Johnson and SAAO-Glass IR-photometry was transformed 
to the SAAO-Carter system according to Glass (1985). 

Table |3] lists the adopted metallicity value. They were 
taken preferentially from Luck & Andrievsky (2004), Kovtyukh 
et al. (2005), and Andrievsky et al. (2005) who studied 
the fundamental parameters, including abundances, at several 
phases. Alternatively, values were taken from other papers by 
Andrievsky /Luck/Kovtyukh, and averaged when multiple deter- 
minations were available. Thirdly, iron abundances from Fry & 
Carney (1997) and Yong et al. (2006) were transformed to the 
"Andrievsky system" following Laney & Caldwell (2007) and 
Luck et al. (2006). Based on the 12 stars in common between 
Andrievsky et al. (2003) and Mottini et al. (2008) an offset AND 
- MOT = H-0.06 ±0.11 was determined, and based on 6 stars 
in common with Lemasle et al. (2007) an offset AND - LEM = 
H-0.09 ± 0.07. These offsets were applied to the Mottini et al. and 
Lemasle et al. iron abundances to bring them to the "Andrievsky 
system". 



3. The model 

The model is outlined in G07 but will be briefly repeated here. 
The V-, K- and RV data with error bars are fitted with a function 
of the form: 

F{t) = -Fo + 5Z (^i '^^^(^TT t e'f) + Bi cos(27r t e'^)) (1) 
j=i 

where P = e^-^ is the period (in days). Typically, P is deter- 
mined from the fit to the available optical photometry as this 
dataset is usually most extensive. The period is then fixed when 
fitting Eq. 1 to the iiT-band and RV data. 

The determination of the parameters is done using the 
MRQMIN routine (using the Levenberg-Marquardt method) from 
Press et al. (1992) written in Fortran77, which minimises 



= ^(F,;-P(t,))V(^F.)^ 



(2) 



http://crocus.physics.mcmaster.ca/Cepheid/ 

Two stars are no longer considered in the Baade-Wesselink analy- 
ses below: BQ Ser and V367 Set which are classified as double-mode 
Cepheids. 



with Fi the measurement at time t[ which has an error bar ctf; 
Also the reduced x^ is defined: 



,2 _ X 



{n~p) 



(3) 
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Table 1. Sources of V-, K-bmd and RV data. 



Name 


V 


K 


RV 


Data not considered^ 


AQPup 


1,2,3,4,5,84 


6,7,8 


1,10,11,44 


3,-,- JD < 34000, JD > 49700 


BBSgr 


2,3,16,63 


6,24 


20,33,39,44 


-,-,44 


beta Dor 


3,11,16 


6,18 


11,19,20,21,22 


-,-,20 


BFOpli 


2,3,16,23,63 


6,18,24 


10,20,39,44,76 


JD < 36000 


BGLac 


2,3,12,13 


12 


14,15,44 




BNPup 


1,3,5,16 


6,7 


1,17 




CDCyg 


2,3,64 


24 


33,44 




CVMon 


2,3,13,16,23 


6,24 


25,26,44 




delta Cep 


2,3,11,12,13,27 


12,98 


4,9,14,25,28,29 


-,-,9 JD < 43500 


DTCyg 


2,3,4,27,30,31 


24,32 


4,9,28,29,33,34 


-,32,9-^28 JD < 34000 


eta Aql 


2,3,12,13,16,27,35 


12,24 


9,14,20,25,28,29,36,37 


-,-,9 


FFAql 


2,3,11,16,27,30,31 


18,24,32 


9,20,28,29,33,38 


3,32,9 JD < 44000 


FMAql 


2,3,12,13,16,23 


12,24 


14,33,39,44 


JD < 30000 


FN Aql 


2,3,12,13,16,23 


12,24 


14,33,39 


-,-,39 


KNCen 


1,3,5,11,16,40,84 


6 


1,17,41 


-,-,41 JD < 40000, JD > 47000 


KQSco 


1,3,11 


6 


1,42 




ICar 


3,5,11,16,35 


6 


11,19,20,23,43 


5,-,23 JD < 46500 


LSPup 


1,3,11 


6 


1 




RSPup 


2,3,5,11,16,84 


6,8,24 


11,19,25,44 


-,-,44 JD < 42000, JD > 50000 


RTAur 


2,12,27,31,45 


12,32 


9,28,29,33,37,46 


-,-,9 


RUSct 


2,3,16 


6 


26,33,105 




RYSco 


1,2,3,5,11,16 


6 


1,20,44 




RYVel 


1,3,5,11,16 


6 


1,11,17 


JD < 44000, JD > 50000 


RZVel 


1,3,5,11,16 


6 


1,11,20,47 




SMus 


3,11,23,40 


6,24 


10,11,20,22,41,47,48,49,50 


-,-,41 JD < 30000 


SNor 


3,11,23,40 


6,24 


4,10,20,41,47,49,51,52,53 


-,-,41 JD < 30000 


SSge 


2,3,11,12,13,23,27,30,40 


12,24,32 


9,14,28,29,33,37,54,55,56 


-,32,54 


SSSct 


2,3,16,23,63 


24 


39,42,44,63 


JD < 30000 


SUCas 


2,3,9,27,31 


12,32 


4,9,25,28,29,33,57,58,59,60,61,62 


-,32,9 JD < 43000 


SUCyg 


2,3,27,30,31 


24 


9,28,33,106,107,108,109 


-,-,9 JD < 35000 


S Vul 


3,30 


6 


33, 110 


JD < 45400 


SVVul 


2,3,4,9,27,64 


6,12,24 


4,9,14,15,25,28,29,33,39,65 


JD < 45500, JD > 48600 


SWVel 


1,3,5,11,16,84 


6 


1,11 




SXVel 


3,23,40 


6 


10 




SZAql 


2,3,11,12,16,23,64 


6,12,24 


11,14,23 




SZTau 


2,3,4,12,27,31 


6,12,32 


4,9,28,29,33 


-,32,9 JD < 43000, JD > 48000 


T Mon 


1,2,3,11,16,35,64 


6,24,32 


4,9,25,29,33,37,48,65,66,67,68,69,70,7 1 ,72,73 


-,32,70 JD < 44500 


TTAql 


2,3,11,12,16,27,64,74 


12,24 


9,11,14,15,22,28,33,37,39,74 


-,-,74 JD < 44000 


TVel 


3,16,75 


6,18 


10,20,75 


-,18,- 


TVul 


2,3,12,27,31 


12,24,32,98 


4,9,14,28,29 


-,32,9 JD < 43600 


UCar 


1,3,5,11,16,35 


6,24 


1,11,20 




UNor 


1,5,11,16 


6 


1,11,17,41 




USgr 


2,3,4,16,23,35,63 


6,24,32 


4,9,10,11,20,25,28,33,39,47,51,76,77,78 


-,32,9 JD < 37000 


UUMus 


1,5,11,40,84 


6 


1,11 




UVul 


2,3,12,13,23,27 


12 


4,14,33,37,44,79,80 


-,-,44 



and the quantity 

BIC = x^ + (p+l) ln(n), 



(4) 



where p = 2 A^ + 2 is the number of free parameters {p ~2N +1 
when fitting the RV and K light curve). As the number N of 
harmonics to be fitted to the data is a priori not known, one 
could obtain ever better fits (lower x^) by increasing N . The 
Bayesian information criterion (Schwarz 1978) is a formalism 
that penalises this, and N (for the V , K and RV curve indepen- 
dently) is chosen such that BIC reaches a minimum. The number 
of harmonics used varies between 3 and 10 in the optical, 1 and 
5 in the NIR, and 2 and 9 for the RV curves. 

Given the analytical form of Eq. 1, the radial velocity curve 
can be integrated exactly to obtain the variation in radius as a 
function of time (phase): 



Ai?(t, 50) 



t+pse 



(wR - 7) dt, 



where 7 is the systemic velocity, wr the radial velocity, p the 
projection factor and 59 allows for a phase shift between the RV 
curve and the angular diameter variations determined via the SB 
relation. 

Then, the equation 



e{t) =9. 3038 mas 



i?o + Ai?(i, 56) 



(5) 



is fitted with the angular diameter in mas, Rq the stars radius in 
solar radii and d the distance in pc. The fitting is done using the 
linear bi-sector (using the code SIXLIN from Isobe et al. 1990) 
as used and preferred by e.g. Storm et al. (2004), Barnes et al. 
(2005) and Gieren et al. (2005). In some cases a phase range 
around 0.85-0.95 was excluded from the fit. 

In G07 the p-factor was determined chiefly using six 
Cepheids with interferometrically measured angular diameter 
variations and HST based accurate parallaxes (Benedict et al. 
2007). It was concluded that, formally, there was no need for 
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Table 1. Continued 



Name 


V 


K 


RV 


Data not considered^ 


V340 Nor 


3,4,82 


6 


4,26 




V350 Sgr 


2,3,23,63,99 


18,24 


20,22,33,44,76,83,39 


-,24,- JD < 42000 


V496 Aql 


2,3,16,63 


24 


9,10,20,33,76 




VCar 


3,16,23 


6 


10.20,47 




VCen 


3,11,16,23,63 


6,24 


10,11,20,47,76 


-,24,- JD < 42000 


VWCen 


1,3,16,84 


6 


1,41 


-,-,41 


VYCar 


1,3,5,11,16,99,100 


6,24 


1,11,20,41 


JD < 42000 


VZCyg 


2,3,12,31 


12,24 


4,14,33,39,44 


JD < 40000 


VZPup 


1,3,5,11,16 


6,7,8 


1,11 




WSgr 


2,3,40,85 


24,32,86 


4,12,20,22,28,47,48,85,87,88,89 




WZCar 


1,3,5,11,16 


6 


1,11 




WZSgr 


1,2,3,5,11,40 


6 


1,11,33,39 


JD < 40000 


XCyg 


2,3,27 


12,24,32,98 


4,14,25,29,33 




XLac 


2,3,12,13 


12 


4,14,33,44 


-,-,44 


XPup 


2,3,5,11,16,23 


6 


11,23,44 


5,-,- 


XSgr 


2,3,11,16,35,99 


24,32,104 


9,10,11,20,28,47,52,101,102,103 


-,-,9 JD < 35000 


XXCen 


3,11,40,74 


6,24 


10,11,20,41,52,74 


-,24,74 


YLac 


2,3,12,31,90 


12 


14,15,33,39,44 


JD < 26500 


YOph 


1,2,3,16 


6,24 


1,9,19,20,28,33,34,91,92 


-,-,9 JD < 40000 


YSgr 


2,3,11,16,23,35 


24,32 


9,11,19,20,28,46,49,93,94 


-,-,9 JD < 28500 


YZSgr 


2,3,23,40 


24 


20,44 




ZLac 


2,3,12,64 


12 


14,33,39,44,80,95,96,97 


JD < 40500 


zeta Gem 


2,3,4,27,35,64 


32,104 


4,9,19,28,29,33,37 


-,-,9 



References: 

1= Coulson & Caldwell (1985a); 2= Moffett & Barnes (1984); 3= Berdnikov et al. (2000), a datafile named "cepheids- 16-03-2006" was retrieved 
from the ftp address Hsted in that paper; 4= Bersier et al. (1994); 5= Madore (1975); 6= Laney & Stobie (1992); 7= Schechter et al. (1992); 8= 
Welch (1985); 9= Barnes et al. (1987), points with uncertainty flag ":" were removed; 10= Stibbs (1955); 11= Bersier (2002), data points with 
weight and 1 in the Geneva photometry were removed; 12= Barnes et al. (1997); 13= Szabados (1980); 14= Barnes et al. (2005); 15= Imbert 
(1999); 16= Pel (1976); 17= Pont et al. (1994); 18= Lloyd Evans (1980a); 19= Nardetto et al. (2006); 20= Lloyd Evans (1980b); 21= Taylor & 
Booth (1998); 22= Petterson et al. (2005); 23= Caldwell et al. (2001); 24= Welch et al. (1984); 25= Storm et al. (2004); 26= Metzger et al. (1992); 
27= Kiss (1998); 28= Wilson et al. (1989); 29= Kiss (2000); 30= Szabados (1991); 31= Szabados (1977); 32= Wisniewski & Johnson (1968); 33= 
Gorynya et al. (1998, VizieR On-line Data Catalog: III/229); 34= Sanford (1935); 35= Shobbrook (1992); 36= Jacobsen &l Wallerstein (1981); 
37= Evans (1976); 38= Evans et al. (1990); 39= Joy (1937); 40= Walraven et al. (1964); 41= Grayzeck (1978); 42= This paper, TableH 43= 
Taylor et al. (1997); 44= Barnes et al. (1988), points with uncertainty flag ":" were removed; 45= Turner et al. (2007); 46= Duncan (1908); 47= 
Lloyd Evans (1968); 48= Petterson et al. (2004); 49= Campbefl & Moore (1928); 50= Bohm-Vitense et al. (1990); 51= Mermifliod et al. (1987); 
52= Feast (1967); 53= Breger (1970); 54= Evans et al. (1993); 55= Herbig & Moore (1952); 56= Breitfellner & Gillet (1993); 57= Adams & 
Shapley (1918); 58= Abt (1959); 59= Niva & Schmidt (1979); 60= Gieren (1976); 61= Haupl (1988); 62= Beavers & Eitter (1986); 63= Gieren 
(1981b); 64= Szabados (1981); 65= Sanford (1956); 66= Evans et al. (1999); 67= Frost (1906); 68= Sanford (1927); 69= Wallerstein (1972); 70= 
Coulson (1983); 71= Evans & Lyons (1994); 72= Gieren (1989); 73= Harper (1934); 74= Coulson et al. (1985); 75= Gieren (1985); 76= Gieren 
(1981a); 77= Jacobsen (1970); 78= Breger (1967); 79= Sanford (1951); 80= Imbert (1996); 81= Coulson & Caldwell (1985b); 82= Eggen (1983); 
83= Evans & Sugars (1997); 84= Hipparcos Epoch Photometry; 85= Babel et al. (1989); 86= Kimeswenger et al. (2004); 87= Albrow & Cottrell 
(1996); 88= Jacobsen et al. (1984); 89= Jacobsen (1974); 90= Henden (1980); 91= Evans & Lyons (1992); 92= Abt & Levy (1978); 93= Duncan 
(1922); 94= ten Bruggencate (1930); 95= Evans & Welch (1993); 96= Gieren (1989); 97= Sugars & Evans (1996); 98= Fernley et al. (1989); 99= 
Dean (1977); 100= Dean (1981); 101= Duncan (1932); 102= Slipher (1904); 103= Moore (1909); 104= Feast et al. (2008); 105= Metzger et al. 
(1993); 106= Evans (1988); 107= Maddrifl (1906); 108= Abt (1973); 109= Imbert (1985); 110= Joy (1952). 
Note: 1= The number indicates the dataset not used in, respectively, V , K, RV. 



a period-dependent p-factor, and p = 1.27 ± 0.05 fitted the 
available data. On the other hand, a moderate dependence of 
p ~ — O.OSlogP as used in the Gieren et al. papers (1993, 
1997, 1998, Storm et al. 2004, Barnes et al. 2003), or p - 
—0.064 log P as advocated by Nardetto et al. (2007) are also 
consistent with the available data. 

The Nardetto et al. (2007) theoretical investigation suggest 
that there is a difference between the p-factor to be used with 
wide-band interferometry (like in G07) and with RV data (when 
applying the SB technique as in the present study). For 6 Cep this 
difference is of the order of 0.06 (Nardetto et al. 2004, 2007). 

In the present paper, two p-factors will be considered: p = 
1.33 (the value found in G07 corrected by 0.06), and p — 1.376— 
0.064 log P from Nardetto et al. (2007). The derived distances 
depend linearly on the adopted value of p. 



An SB relation can be defined as follows (see van Belle 
1999): 



9o^ex io("i/5), 



(6) 



where 9 is the LD angular diameter (in mas), and mi a de- 
reddened magnitude (for example, Vq). The logarithm of this 
quantity (the zero magnitude angular diameter) is plotted against 
a de-reddened colour (for example, {V — K)o), 



log^o = a X (?7i2 — m:i) + b. 



(7) 



G07 analysed the available datasets of Cepheids with inter- 
ferometrically determined angular diameters over the light curve 
and derived a SB relation. In a recent paper. Feast et al. (2008) 
presented new NIR data for zeta Gem and X Sgr Taking this 
new data into account, Eq. 9 in G07 becomes 

log 00 = (0.2724 ± 0.0045) {V - K)o + (0.5283 ± 0.0091)(8) 
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Table 2. New radial velocity data 



Table 3. Adopted [Fe/H] values for the Cepheids in the sample 



JD 



RV (km s~') 



JD 



RV(kms"') 



KQSco 




V367 Set 




2454302.678 


-7.57 


2454341.605 


-22.61 


2454303.670 


-10.77 


2454341.605 


26.97 


2454304.527 


-26.89 


2454302.737 


-19.57 


2454307.627 


-45.12 


2454303.709 


-19.15 


2454308.672 


-45.51 


2454305.682 


0.77 


2454309.722 


-53.38 


2454306.754 


4.38 


2454310.570 


-49.50 


2454307.780 


-16.31 


2454311.711 


-48.46 


2454308.741 


-16.90 


2454312.605 


-50.97 


2454310.747 


-5.25 


2454332.497 


-12.06 


2454311.733 


-7.84 


2454336.489 


-47.39 


2454312.762 


-9.21 


2454339.524 


-42.66 


2454332.588 


5.21 


2454341.537 


-46.58 


2454334.567 


-20.65 






2454334.581 


-20.30 


SSSct 




2454339.626 


-13.20 


2454302.580 


7.55 


2454341.580 


-9.54 


2454303.739 


-24.37 






2454304.754 


-11.28 


BQSer 




2454305.772 


3.70 


2454302.601 


-18.97 


2454306.711 


1.72 


2454303.731 


-11.41 


2454307.691 


-23.22 


2454304.724 


-15.08 


2454308.624 


-8.05 


2454305.738 


-11.01 


2454309.616 


5.48 


2454306.731 


-1.57 


2454310.785 


-17.19 


2454307.747 


-27.11 


2454311.656 


-18.74 


2454308.717 


-17.00 


2454312.806 


-0.36 


2454309.737 


-4.14 


2454331.699 


6.05 


2454310.705 


-2.98 


2454332.556 


-4.74 


2454311.753 


-22.50 


2454334.613 


-3.38 


2454312.718 


-11.03 


2454335.533 


7.01 


2454331.653 


1.03 


2454336.671 


-23.15 


2454332.517 


-23.78 


2454339.644 


5.00 


2454334.536 


-10.36 


2454341.637 


-8.36 


2454335.515 


-13.03 






2454339.557 


-6.35 






2454341.559 


-25.73 



which is the form of the the SB relation used in this paper 

Colour excess values E{B — V) used to de-redden the V and 
K magnitudes have been taken from the compilation by Fouque 
et al. (2007; hereafter F07). Reddening ratios of Ry= 3.3 and 
^k/^v= 0.0909 have been adopted, as in G07. For the stars 
already studied in G07 (beta Dor, 1 Car, zeta Gem, W Sgr, Y Sgr, 
del Cep, FF Aql, T Vul, RT Aur), for consistency, the E{B - V) 
values from G07 have been taken, which are based on the Ay 
values in Benedict et al. (2007) and the reddening ratio of 3.3. 
The resulting E{B — V) values do not differ significantly from 
the values listed in F07. 



4. Binary Cepheids 

A number of stars in the sample are known or suspected spec- 
troscopic binaries. In order to apply the SB technique outlined 
in the previous section, the RV data have to be corrected for the 
binary motion. 

The database on Cepheids in binary system^ (Szabados 
2003) is most useful as it lists published orbital parameters and 
references. For some stars, more RV data have become available 
than was used in the last paper that published an orbit and so it 
was decided to rederive the 5 orbital parameters (in the case of 
the eccentricity. In e was the parameter actually fitted). 



Name 


[Fe/H] 


Name 


[ Fe/H] 


AQPup 


-0.14(5) 


SZAql 


-1-0.17 ±0.04 (2) 


BBSgr 


-^0.08 (6) 


SZTau 


-1-0.08 (4) 


beta Dor 


-0.01 (4) 


TMon 


-1-0.15 ±0.03 (2) 


BFOph 


-hO.OO (6) 


TTAql 


-1-0.10 ±0.02 (2) 


BGLac 


-0.01 (4) 


TVel 


-0.02 (5) 


BNPup 


-^0.01 (4) 


TVul 


-1-0.01 ± 0.02 (3) 


CDCyg 


-^0.11 ±0.05(2) 


UCar 


-1-0.23 (7) 


CVMon 


-0.03 (4) 


UNor 


-1-0.13(7) 


del Cep 


-1-0.08 (4,6) 


USgr 


-1-0.08 ± 0.02 (1) 


DTCyg 


-^0.11 (4) 


UUMus 


-1-0.17(7) 


eta Aql 


-^0.08 ± 0.02 (1) 


UVul 


-1-0.11 ±0.03(1) 


FFAql 


-1-0.02 (4) 


V340 Nor 


-1-0.00 (4) 


FMAql 


-fO.08 (4) 


V350 Sgr 


-1-0.18(6) 


FN Aql 


-0.02 (4) 


V496 Aql 


-1-0.05 (4) 


KNCen 


-hO.23 (7) 


VCar 


-1-0.02 (7,8) 


KQSco 


-h0.16(4) 


VCen 


-1-0.04 (4) 


ICar 


-^0.13 (7,8) 


VWCen 


-1-0.04 (7) 


LSPup 


-0.10(7) 


VYCar 


-1-0.03 (8) 


RSPup 


-^0.17 (5) 


VZCyg 


-1-0.05 (6) 


RTAur 


-1-0.06 (4) 


VZPup 


-0.15(4) 


RUSct 


-1-0.01 (4) 


WSgr 


-1-0.03 ± 0.03 (1) 


RYSco 


-hO.09 (6) 


WZCar 


-1-0.14(7) 


RYVel 


-0.03 (5) 


WZSgr 


-1-0.20 ± 0.05 (2) 


RZVel 


-0.07 (5) 


XCyg 


-1-0.10 ±0.03 (2) 


SMus 


-1-0.19(7) 


XLac 


-0.02 (6) 


SNor 


-^0.05 (4) 


XPup 


-0.03 ± 0.03 (2) 


SSge 


-1-0.10 ±0.02(1) 


XSgr 


-0.29 (4) 


SSSct 


-1-0.06 (6) 


XXCen 


-1-0.10(7) 


SUCas 


-0.01 (4) 


YLac 


-0.04 ± 0.06 (3) 


SUCyg 


-0.03 ± 0.04 (3) 


YOph 


-1-0.05 (4) 


SVul 


-0.03 ± 0.03 (2) 


YSgr 


-1-0.05 ± 0.03 (3) 


SVVul 


-^0.05 ± 0.04 (2) 


YZSgr 


-1-0.08 (1) 


SWVel 


-0.03 (4,5) 


zeta Gem 


-1-0.04 (4) 


SXVel 


-0.03 (5) 


ZLac 


-1-0.01 ± 0.03 (2) 



http://www.konkoly.hu/CEP/intro.html 



1= Luck & Andrievsky (2004); 2= Kovtyukh et al. (2005b); 3= 
Andrievsky et al. (2005); 4= Andrievsky et al. (2003); 5= Luck et al 
(2003); 6= Luck et al. (2006); 7= Mottini et al. (2008); 8= Lemasle et 
al. (2007). 



A weighted fit to the RV curve was performed using 
MRQMIN with (2 N H- 2) H- 5 free parameters (sometimes less 
when e.g. the eccentricity was fixed), with N the number of 
Fourier terms in the pulsation curve. This procedure does not 
take into account the possibility that the period also changed 
over the time spanned by the RV data. 

As the data sometimes spans more than a hundred years of 
data with possible shifts in zero points and varying intrinsic ac- 
curacies, the fitting was done as follows. Initially all data points 
were given the same error in RV, and the fit was made. A mean 
offset and rms per dataset between data and model was calcu- 
lated. If the offset was more than half the rms, the offset was ap- 
plied. A second fit was made. The binary motion was subtracted 
and the datasets compared as outlined above in order to assign 
realistic error bars per dataset. The fit to the binary + pulsation 
RV curve was repeated to arrive at the final fit. 

The orbital elements are listed in Table |4] For most of the 
known spectroscopic binaries the derived orbital parameters are 
close to literature values, cf. Evans et al. (1990) for FF Aql, 
Imbert (1996) for Z Lac and U Vul, Evans et al. (1999) for T 
Mon, Petterson et al. (2004) for S Mus and V350 Sgr, Benedict 
et al. (2007) for W Sgr, Gorynya et al. (1995) for S Sge, Evans 
(1988) for SUCyg. 
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The only exceptions are VZ Cyg and SU Cas. For VZ Cyg, 
Gorynya et al. (1995) quote an orbital period of 725 days, a 
velocity amplitude of 6.5 km s^^ and an eccentricity of 0.05, 
while Rastorgouev et al. (1997) give very different values of 
P= 1483d, K= 15.1 km s~\ e= 0.74. After pre-whitening the 
RV curve with the pulsation period, a Fourier analysis using 
Period04 indicates a period of 2080 days with an amplitude of 
3.3 km s^^. Using a larger dataset than the previous papers the 
eccentric orbit could not be confirmed, and Table |4] lists the de- 
rived elements for e= 0, which results in an orbital period of 
2092 days. For SU Cas, Gorynya et al. (1996) quote a period of 
408 days, velocity amplitude of 3 km s^^ and e= 0.43, while 
Rastorgouev et al. (1997) give P= 407d, K= 5.2 km s~\ e= 
0.73. After pre-whitening the RV curve with the pulsation pe- 
riod, a Fourier analysis using Period04 indicates a period of 408 
days with an amplitude of 1.2 km s^^. Using the larger dataset 
than the previous works, the eccentric orbit could not be con- 
firmed, and Table |4] lists the derived elements for e= 0, which 
results in an orbital period of 407 days and a small velocity am- 
plitude. 

For four stars no orbital elements seem to have been pub- 
lished before. For V496 Aql, Szabados (1989) suggested peri- 
ods of 1200, 1780, 2700, 3600, 5370, or 10750 days. After pre- 
whitening the RV curve with the pulsation period a Fourier anal- 
ysis using Period04 indicates periods of 529, 1355, 1825, 4065 
days with an amplitude of 2.6-3.4 km s^^. 

For a fixed velocity amplitude of 3 km s~^ and zero eccen- 
tricity the period parameter space was searched for the best fit. 
This was found for a period of 1331 days, closely followed by 
a fit with a period of 1769 days. Other periods resulted in much 
poorer fits. Although clearly more RV data are needed to bet- 
ter determine the orbital parameters, the RV curve was corrected 
using the parameters listed in Table 2] 

For X Sgr a period significantly longer than the 507.25 days 
suggested by Szabados (1990) is derived even though only a few 
additional RV have become available since. The orbital elements 
in Table|4]assume zero eccentricity. 

For S Nor, Szabados (1989) suggested periods of 3300 or 
6500 days. After pre-whitening the RV curve with the pulsation 
period a Fourier analysis using Period04 indicates a period of 
7910 days with an amplitude of 3.2 km s^^. Fixing the eccen- 
tricity at zero, the best fitting period is found to be near 3600 
days. Although more RV data are needed to better determine the 
orbital parameters, the RV curve was corrected using the param- 
eters listed in Table |4] 

For XX Cen, Szabados (1990) suggested a period of 909 
days which is largely confirmed. Although clearly more RV data 
are needed to better determine the orbital parameters, the RV 
curve was corrected using the parameters listed in Table|4]which 
have been derived for zero eccentricity. 

Finally, for Y Oph the suggestion in Evans & Lyons (1986) 
that there is little evidence for orbital motion is confirmed us- 
ing significantly more data than they had at their disposal. For 
LS Pup, Szabados (1996) suggested it to be a binary but no new 
RV data have been published for this star in the last ten years 
and no correction for binary motion has been attempted. For Y 
Sgr, Szabados (1989) suggested a binary period of order 10000- 
12000 days (and longer could not be excluded). Assuming zero 
excentricity, several oribital periods longer than 10000 days 
were tried, but no convincing solution could be found. For BF 
Oph, Szabados (1989) suggested a binary period of 4420 ± 80 
days. After pre-whitening the RV curve with the pulsation pe- 
riod, a Fourier analysis using Period04 indicates a period of 4085 
days with a velocity amplitude near 3 km s^^. The least-squares 



fitting to the RV data then revealed an alternative good fit with a 
period of 24700 ± 2900 days with an amplitude of 2.5 km s~^. 
In the end the RV data of BF Oph was not corrected. 



5. Results 

5.1. Distances and magnitudes 

Table |5] lists the distances, radii and absolute magnitudes ob- 
tained. The table also lists the adopted E{B — V) and error 
bar, the derived period, and the p-factor following p — 1.376 — 
0.064 log P from Nardetto et al. (2007). The distances, radii and 
errors scale directly with p, so the values for p = 1.33 are not 
repeated. The error in the period is a few units in the last deci- 
mal place. For the derived quantities two error bars are quoted. 
For the distance and radius the first error bar listed is the error in 
the fit, and for the absolute magnitudes the error due to the er- 
ror in distance and E{B — V). The second error quoted is based 
on a Monte Carlo simulation where new datasets are generated 
based on the error bar in each individual V, K, RV measure- 
ment, and analysed taking an E{B — V) value randomly drawn 
from a Gaussian distribution based on the listed mean value and 
la error bar The second error quoted is the Icr dispersion in the 
derived quantities. 

Figure [T] illustrates the fit to the V, K, and RV curve for AQ 
Pup, while Fig. |2] shows the variation of the angular diameter 
against phase and the change in angular diameter derived from 
the SB relation against the change in radius from integration of 
the RV curve from which the distance is derived (see Eq. |5]l. 
Figures similar to[T]and|2]for all stars in the sample are available 
from the author. 

For ten stars the best available distance comes from the HST- 
based parallax (Benedict et al. 2007). Although a BW-type anal- 
ysis has been carried out for these stars (also see G07), only the 
mean V and K magnitudes derived here have been used together 
with the parallax from Benedict et al. to derive the absolute mag- 
nitudes listed in Table |6] 

F07 studied the calibration of the Galactic PL-relation in 
various photometric bands using a set of 59 calibrating stars, 
with distances from HST and revised Hipparcos parallaxes, in- 
terferometric Baade-Wesselink and infrared SB parallaxes and 
Main-Sequence fitting. 

Of the 58 stars in Table|5] 57 have a distance quoted in F07. 
The ratio of the distance quoted in Table |5] and that in Table 6 
of finally adopted parallaxes in F07 is 1 .02 ± 0.11 (minus one 
outlier, X Lac). The absolute difference in distance between the 
two estimates is 1 . 1 times the combined error bar. Regarding the 
absolute magnitudes quoted in their Table 7, one difference w.r.t. 
the values in Table |5] is due to the fact that F07 use published 
intensity-mean values, while in the present paper these have been 
derived in a self-consistent way from the Fourier fit to the data. 

5.2. PL{Z) relations 

Figures [3] and |4] show the Pi-relation in the F-band and K- 
band assuming the Nardetto et al. calibration of the p-factor, and 
in the bottom panel the residual plotted against metallicity. The 
first and second error bar quoted for the magnitudes in Table |5] 
have been added in quadrature. The period is the fundamental 
one, and the classification into fundamental and first overtone 
pulsator is taken from F07. Four stars in the sample are overtone 
pulsators (DT Cyg, SU Cas, SZ Tau, X Lac) and their period is 
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Table 4. Derived orbital parameters of binary Cepheids. Quantities without error bar have been fixed. 



Name 



K(kms-') 



0.027 ±-t^ 
0.086 i«;rr 



n 9^8 +0.005 
U./J8 -0.005 



^ (degr) 



To (JD-2400000) Period (d) 



FFAql 

SMus 

SNor 

SSge 

SUCas 

SUCyg 

TMon 

UVul 

V350 Sgr 

V496 Aql 

VZCyg 

WSgr 

XSgr 

XX Can 

Z Lac 



4.91 ± 0.07 
14.78 ±0.10 
2.53 ± 0.35 
15.62 ± 0.06 
0.98 ± 0.08 
29.83 ±0.15 

8.42 ±0.19 
3.64 ± 0.44 
11.21 ±0.12 

3.0 

3.46 ±0.19 

1.43 ± 0.06 
2.30 ± 0.27 

4.47 ± 0.28 
10.44 ±0.10 







0.350 
0.414 
0.675 
0.405 



+0.004 
-0.003 
+0.013 
-0.013 
+0.033 
-0.031 
+0.015 
0.014 








0.153 



+0.046 






O.UZJ _o,oo8 



319 ±45 


45437 ± 178 


1432.4 ± 1.1 


194 ±4 


48575 ± 5 


505.15 ±0.08 





45639 ± 67 


3584 ± 33 


202.7 ±1.1 


48010.7 ± 1.9 


675.72 ± 0.04 





50279 ± 6 


406.76 ± 0.04 


223.85 ± 0.85 


43766.2 ± 1.1 


549.24 ± 0.02 


204 ±3 


49300 ± 143 


32449 ± 726 


353 ±4 


44800 ± 16 


2510 ±2.8 


279 ±2 


47594 ± 9 


1482 ± 2.4 





45606 ± 25 


1331 ±6.5 





45002 ± 50 


2092 ± 21 


331 ± 16 


48380 ± 79 


1651 ±8.8 





48208 ± 19 


573.6 ± 0.6 





44861 ± 8 


924.1 ± 1.1 


344 ± 21 


46582 ± 22 


382.63 ±0.10 




Fig. 1. The phased curves in V, K and RV are shown for AQ 
Pup. Data points are shown with errors bars and the line shows 
the harmonic fit. The fits to all stars in the sample are available 
from the author upon request. 




-10 

AR (solar radii) 




Phase 

Fig. 2. For AQ Pup, the top panel shows the linear-bisector fit 
to the angular diameter as a function of radial displacement. The 
bottom panel shows the angular diameter against phase. Crosses 
represent datapoints not considered in the fit. The fits to all stars 
in the sample are available from the author upon request. 



fundamentalised using an iterative procedure using the formula 
(Alcocketal. 1995) 



Pi/Pq ^ 0.720 - 0.027 log Po 



(9) 



The metallicity dependence has been determined in two 
ways, one by first fitting a linear PL-relation, and fitting the 
residual with a linear relation against [Fe/H] (as shown in Figs.|3] 
andlU, and secondly by making a linear-fit in the two variables 
logP and [Fe/H], as quoted below. The results are essentially 
the same. 
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Table 5. Distances, radii and absolute magnitudes from the BW-analysis 



Name 


E{B - V) 


Period (d) 


P 


d 


-R(Rq) 


My 


Mk 


AQPup 


0.518 ±0.010 


30.0980 


1.281 


3292.2 ± 80.0 ± 125.8 


149.7 ± 3.7 ± 6.0 


-5.51 ±0.06 ±0.09 


-7.44 ± 0.05 ± 0.09 


BBSgr 


0.281 ±0.009 


6.637111 


1.323 


768.7 ± 25.0 ± 25.7 


46.9 ± 1.5 ± 1.5 


-3.40 ± 0.08 ± 0.08 


-5.00 ± 0.07 ± 0.07 


BFOph 


0.235 ± 0.010 


4.067677 


1.337 


727.7 ± 32.6 ± 28.9 


31.9 ± 1.4 ± 1.3 


-2.72 ±0.10 ±0.09 


-4.20 ±0.10 ±0.09 


BGLac 


0.300 ± 0.016 


5.33192 


1.329 


1700.9 ± 24.3 ± 38.5 


39.9 ± 0.6 ± 0.9 


-3.24 ± 0.06 ± 0.08 


-4.70 ± 0.03 ± 0.05 


BNPup 


0.416 ±0.018 


13.67243 


1.303 


3773.4 ± 94.6 ± 145.9 


80.2 ± 2.0 ± 3.2 


-4.30 ±0.08 ±0.11 


-6.08 ± 0.05 ± 0.09 


CDCyg 


0.493 ± 0.015 


17.07406 


1.297 


2336.7 ±64.1 ±60.5 


88.3 ± 2.4 ± 2.2 


-4.46 ± 0.08 ± 0.08 


-6.31 ±0.06 ±0.06 


CVMon 


0.722 ± 0.022 


5.378662 


1.329 


1475.1 ±60.7 ±53.9 


36.7 ± 1.5 ± 1.3 


-2.92 ±0.11 ±0.12 


-4.48 ± 0.09 ± 0.08 


DTCyg 


0.042 ±0.011 


2.499189 


1.351 


603.3 ± 75.7 ± 66.5 


34.9 ± 4.4 ± 3.9 


-3.26 ± 0.26 ± 0.29 


-4.48 ± 0.26 ± 0.28 


eta Aql 


0.130 ±0.009 


7.176813 


1.321 


261.0 ±6.7 ±7.1 


49.3 ± 1.3 ± 1.3 


-3.59 ± 0.06 ± 0.07 


-5.14 ±0.05 ±0.06 


FMAql 


0.589 ±0.012 


6.114274 


1.326 


1141.2 ±38.8 ±37.4 


55.1 ± 1.9 ± 1.8 


-3.93 ± 0.08 ± 0.09 


-5.40 ± 0.07 ± 0.07 


FN Aql 


0.483 ± 0.010 


9.48163 


1.313 


1217.5 ±28.8 ±25.9 


52.2 ± 1.2 ± 1.1 


-3.62 ± 0.06 ± 0.06 


-5.21 ± 0.05 ± 0.04 


KNCen 


0.797 ± 0.091 


34.0498 


1.278 


3973.5 ± 93.4 ± 161.9 


177.1 ±4.2 ±4.5 


-5.71 ±0.30 ±0.35 


-7.74 ±0.06 ±0.11 


KQSco 


0.869 ± 0.021 


28.6958 


1.283 


2655.2 ± 56.4 ± 105.8 


160.9 ± 3.4 ± 6.2 


-5.13 ±0.08 ±0.12 


-7.43 ± 0.05 ± 0.09 


LSPup 


0.461 ± 0.015 


14.14725 


1.302 


4919.8 ± 91.3 ± 157.8 


85.4 ± 1.6 ±2.7 


-4.49 ± 0.06 ± 0.09 


-6.24 ± 0.04 ± 0.07 


RSPup 


0.457 ± 0.009 


41.457 


1.272 


1765.4 ± 33.9 ± 65.2 


174.2 ±3.4 ±6.8 


-5.67 ± 0.05 ± 0.09 


-7.71 ± 0.04 ± 0.08 


RUSct 


0.921 ±0.012 


19.7032 


1.293 


2043.5 ± 78.3 ± 60.9 


107.6 ±4.1 ±3.2 


-5.04 ± 0.09 ± 0.08 


-6.76 ± 0.08 ± 0.07 


RYSco 


0.718 ±0.018 


20.32134 


1.292 


1185.1 ±26.0 ±56.8 


94.1 ±2.1 ±4.4 


-4.68 ±0.08 ±0.13 


-6.44 ±0.05 ±0.11 


RYVel 


0.547 ± 0.010 


28.1240 


1.283 


2178.5 ±28.1 ±95.3 


118.6 ± 1.5 ±5.0 


-5.08 ±0.04 ±0.10 


-6.92 ± 0.03 ± 0.09 


RZVel 


0.299 ± 0.010 


20.3996 


1.292 


1472.3 ± 17.9 ±27.1 


103.8 ± 1.3 ± 1.8 


-4.67 ± 0.04 ± 0.06 


-6.61 ± 0.03 ± 0.04 


SMus 


0.212 ±0.016 


9.65997 


1.313 


787.7 ± 29.5 ± 24.2 


59.0 ± 2.2 ± 1.8 


-4.03 ±0.10 ±0.09 


-5.53 ± 0.08 ± 0.07 


SNor 


0.179 ±0.009 


9.75425 


1.313 


819.8 ±21.9 ±23.7 


60.0 ± 1.6 ± 1.8 


-3.71 ±0.06 ±0.07 


-5.46 ± 0.06 ± 0.06 


SSge 


0.100 ±0.010 


8.38207 


1.317 


662.9 ± 17.1 ± 14.0 


55.5 ± 1.4 ± 1.2 


-3.80 ± 0.06 ± 0.06 


-5.37 ± 0.06 ± 0.05 


SSSct 


0.325 ± 0.010 


3.671331 


1.340 


959.1 ±90.6 ±71.7 


30.7 ± 2.9 ± 2.4 


-2.76 ±0.20 ±0.18 


-4.15 ±0.20 ±0.17 


SUCas 


0.259 ± 0.010 


1.949329 


1.357 


356.8 ± 18.8 ± 22.9 


23.9 ± 1.3 ± 1.6 


-2.65 ±0.12 ±0.16 


-3.71 ±0.11 ±0.15 


SUCyg 


0.098 ± 0.014 


3.845552 


1.339 


842.1 ±62.0 ±39.0 


32.9 ± 2.4 ± 1.5 


-3.05 ±0.16 ±0.11 


-4.34 ±0.15 ±0.10 


S Vul 


0.727 ± 0.042 


68.714 


1.258 


4430.4 ± 147.7 ± 205.7 


304.8 ± 10.2 ± 12.7 


-6.65 ±0.16 ±0.19 


-8.85 ±0.07 ±0.10 


SVVul 


0.461 ± 0.022 


45.028 


1.270 


2449.9 ± 57.8 ± 58.2 


211.1 ±5.0 ±7.8 


-6.21 ±0.09 ±0.10 


-8.15 ±0.05 ±0.06 


SWVel 


0.344 ± 0.010 


23.4395 


1.288 


2285.6 ± 25.5 ± 54.9 


104.1 ± 1.2 ±2.5 


-4.73 ± 0.04 ± 0.07 


-6.65 ± 0.02 ± 0.05 


SXVel 


0.236 ± 0.012 


9.55042 


1.313 


1788.9 ± 70.4 ± 121.5 


54.5 ±2.1 ±3.7 


-3.74 ±0.09 ±0.15 


-5.33 ±0.08 ±0.15 


SZAql 


0.537 ± 0.017 


17.13977 


1.297 


1960.2 ± 24.3 ± 47.3 


96.4 ± 1.2 ±2.3 


-4.54 ± 0.06 ± 0.09 


-6.46 ± 0.03 ± 0.06 


SZTau 


0.295 ±0.011 


3.148921 


1.344 


525.3 ± 12.1 ±21.6 


33.9 ± 0.8 ± 1.4 


-3.05 ±0.06 ±0.10 


-4.38 ± 0.05 ± 0.09 


T Mon 


0.181 ±0.011 


27.0318 


1.284 


1239.0 ± 20.0 ± 19.6 


125.7 ± 2.0 ± 1.8 


-4.88 ± 0.05 ± 0.06 


-6.98 ± 0.03 ± 0.04 


TTAql 


0.438 ±0.011 


13.75502 


1.303 


1025.1 ±22.1 ±27.2 


82.0 ± 1.8 ±2.1 


-4.32 ± 0.06 ± 0.07 


-6.12 ±0.05 ±0.06 


TVel 


0.289 ± 0.010 


4.639811 


1.333 


942.2 ± 9.2 ± 27.5 


33.9 ±0.3 ±1.0 


-2.78 ± 0.04 ± 0.07 


-4.31 ±0.02 ±0.06 


UCar 


0.265 ± 0.010 


38.8229 


1.274 


1684.2 ±51.9 ±42.8 


170.0 ± 5.3 ± 4.5 


-5.66 ± 0.07 ± 0.07 


-7.68 ± 0.07 ± 0.06 


UNor 


0.862 ± 0.024 


12.64418 


1.305 


1438.3 ± 26.0 ± 32.0 


78.3 ± 1.4 ± 1.6 


-4.37 ±0.09 ±0.10 


-6.05 ± 0.04 ± 0.05 


USgr 


0.403 ± 0.009 


6.745306 


1.323 


541.0 ± 23.9 ± 22.4 


42.9 ± 1.9 ± 1.8 


-3.28 ±0.10 ±0.10 


-4.83 ±0.09 ±0.10 


UUMus 


0.399 ± 0.015 


11.63613 


1.308 


2906.1 ±63.2± 112.0 


64.6 ± 1.4 ±2.5 


-3.79 ±0.07 ±0.10 


-5.60 ± 0.05 ± 0.09 


UVul 


0.603 ±0.011 


7.99074 


1.318 


638.4 ±23.6 ±21.3 


50.6 ± 1.9 ± 1.7 


-3.86 ± 0.09 ± 0.08 


-5.25 ± 0.08 ± 0.07 


V340 Nor 


0.321 ±0.018 


11.2884 


1.309 


1782.4 ± 101.6 ± 131.9 


70.5 ± 4.0 ± 5.2 


-3.90 ±0.13 ±0.18 


-5.76 ±0.12 ±0.17 


V350 Sgr 


0.299 ± 0.010 


5.154248 


1.330 


807.8 ± 18.9 ± 38.0 


35.4 ± 0.8 ± 1.6 


-3.02 ±0.06 ±0.11 


-4.46 ±0.05 ±0.10 


V496 Aql 


0.397 ± 0.010 


6.80702 


1.323 


937.9 ±80.0 ±81.7 


45.5 ± 3.9 ± 4.0 


-3.40 ±0.18 ±0.21 


-4.94 ±0.18 ±0.20 


VCar 


0.169 ±0.011 


6.696707 


1.323 


908.5 ± 28.9 ± 47.0 


38.2 ± 1.2 ±2.0 


-2.96 ±0.08 ±0.12 


-4.56 ±0.07 ±0.12 


VCen 


0.292 ± 0.012 


5.494010 


1.329 


672.6 ± 16.5 ± 20.7 


40.2 ± 1.0 ± 1.2 


-3.25 ± 0.07 ± 0.08 


-4.72 ± 0.05 ± 0.07 


VWCen 


0.428 ± 0.024 


15.03733 


1.301 


3701.4 ± 66.8 ± 120.1 


87.7 ± 1.6 ±2.7 


-3.96 ±0.09 ±0.12 


-6.15 ±0.04 ±0.07 


VYCar 


0.237 ± 0.009 


18.9050 


1.294 


1455.1 ±34.2 ±36.6 


81.6 ± 1.9 ±2.0 


-4.11 ±0.06 ±0.06 


-6.07 ± 0.05 ± 0.06 


VZCyg 


0.266 ±0.011 


4.86440 


1.332 


1773.2 ±51.2 ±54.5 


37.4 ± 1.1 ± 1.2 


-3.13 ±0.07 ±0.08 


-4.57 ± 0.06 ± 0.07 


VZPup 


0.459 ±0.011 


23.17494 


1.289 


4315.0 ± 86.7 ± 162.2 


99.7 ± 2.0 ± 4.0 


-4.98 ±0.06 ±0.10 


-6.63 ± 0.04 ± 0.09 


WZCar 


0.370 ±0.011 


23.01521 


1.289 


3600.9 ± 50.8 ± 55.7 


101.4 ± 1.4 ± 1.8 


-4.66 ± 0.05 ± 0.06 


-6.59 ± 0.03 ± 0.04 


WZSgr 


0.431 ±0.011 


21.8501 


1.290 


1754.1 ±29.8 ±26.8 


117.6 ± 2.0 ± 1.7 


-4.55 ± 0.05 ± 0.06 


-6.77 ± 0.04 ± 0.04 


XCyg 


0.228 ± 0.012 


16.3856 


1.298 


1081.0 ±25.3 ±22.7 


92.8 ± 2.2 ± 1.9 


-4.49 ± 0.06 ± 0.07 


-6.39 ± 0.05 ± 0.05 


XLac 


0.336 ±0.011 


5.44453 


1.329 


1507.3 ± 83.0 ± 70.9 


44.3 ±2.4 ±2.1 


-3.59 ±0.12 ±0.12 


-4.95 ±0.12 ±0.11 


XPup 


0.402 ± 0.009 


25.9657 


1.285 


2779.4 ± 48.0 ± 133.1 


117.0 ±2.0 ±5.8 


-4.95 ±0.05 ±0.11 


-6.90 ±0.04 ±0.10 


XXCen 


0.266 ±0.011 


10.95351 


1.309 


1412.5 ±32.9 ±31.9 


58.0 ± 1.4 ± 1.2 


-3.77 ± 0.06 ± 0.06 


-5.42 ± 0.05 ± 0.05 


YLac 


0.207 ± 0.016 


4.323760 


1.335 


2523.9 ±74.1 ±86.0 


41.8 ± 1.2 ±1.4 


-3.52 ±0.08 ±0.10 


-4.85 ± 0.06 ± 0.07 


YOph 


0.645 ± 0.015 


17.12614 


1.297 


594.4 ± 17.5 ± 21.0 


85.8 ± 2.5 ± 3.0 


-4.84 ±0.08 ±0.10 


-6.33 ± 0.06 ± 0.08 


YZSgr 


0.281 ±0.010 


9.55376 


1.313 


1102.4 ±53.9 ±90.8 


57.2 ±2.8 ±4.8 


-3.77 ±0.11 ±0.19 


-5.41 ±0.10 ±0.19 


ZLac 


0.370 ±0.011 


10.88569 


1.310 


1877.5 ±37.4 ±40.1 


68.4 ± 1.4 ± 1.5 


-4.14 ±0.06 ±0.07 


-5.78 ± 0.04 ± 0.05 



For the adopted p-factor dependence on period the PLZ- and in K 
relation in V becomes 



Mv = (-2.60 ± 0.09) logP ± (-1.30 ± 0.10) 

±(±0.27±0.30)[Fc/H], (10) 



Mk = (-3.38 ± 0.08) logP ± (-2.19 ± 0.09) 

±(-0.11±0.24)[Fc/H]. (11) 
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Table 6. Absolute magnitudes for the Cepheids with an HST based parallax 



Name 


E{B - V) 


TT 


Mv 


Mk 


beta Dor 


0.076 ± 0.015 


3.14 ±0.16 


-4.01 ±0.12 


-5.58 ±0.11 


del Cep 


0.070 ± 0.010 


3.66 ±0.15 


-3.43 ± 0.09 


-4.87 ± 0.09 


FFAql 


0.194 ±0.018 


2.18 ±0.18 


-3.02 ±0.15 


-4.32 ±0.14 


ICar 


0.158 ±0.018 


2.01 ± 0.20 


-5.25 ± 0.22 


-7.43 ± 0.21 


RTAur 


0.061 ± 0.024 


2.40 ±0.19 


-2.82 ±0.18 


-4.19 ±0.17 


TVul 


0.103 ±0.018 


1.90 ±0.23 


-3.18 ±0.26 


-4.43 ± 0.25 


WSgr 


0.112 ±0.009 


2.28 ± 0.20 


-3.90 ±0.19 


-5.43 ±0.18 


XSgr 


0.176 ±0.030 


3.00 ±0.18 


-3.62 ±0.16 


-5.14 ±0.13 


YSgr 


0.203 ± 0.012 


2.13 ±0.29 


-3.26 ± 0.28 


-4.80 ± 0.28 


zeta Gem 


0.018 ±0.010 


2.78 ±0.18 


-3.94 ±0.14 


-5.65 ±0.14 



For a constant p-factor of 1.33 this would have become 
Mv = (-2.72 ± 0.09) logP± (-1.23 ±0.10) 

±(±0.24±0.30)[Fc/H], (12) 
and in K 
Mk = (-3.50 ± 0.08) logP ± (-2.09 ± 0.09) 

±(-0.10±0.24)[Fc/H]. (13) 

Ignoring the metallicity dependence results in PL-relations that 
have slopes and zeropoints that differ in the third decimal from 
the ones quoted above. 

Ignoring RS Pup, SV Vul and S Vul which have the longest 
periods, and hence limiting the period to <40 days results in 

Mv = (-2.50 ± 0.10) logP± (-1.40 ±0.11) 

±(±0.14±0.31)[Fc/H], (14) 
and in K 
Mk = (-3.29 ± 0.08) logP ± (-2.28 ± 0.10) 

±(-0.01±0.26)[Fe/H]. (15) 
For a constant p-factor of 1.33 this would have become 
Mv = (-2.60 ± 0.10) logP± (-1.33 ±0.11) 

±(±0.11±0.31)[Fc/H], (16) 
and in K 
/Mk = (-3.41 ± 0.08) logP ± (-2.18 ± 0.10) 

±(-0.02±0.25)[Fc/H]. (17) 

The effect of the negative log P dependence on p compared 
to a constant value is to flatten the PL-relation. The size of the 
effect is expected to be about 5 log (1 + 0.064/1.33) « 0.10 and 
this is indeed reflected in the fitted slope. This change in slope is 
about equal to the la error bar derived in the slope. 

F07 quotes slopes of -2.77 ± 0.08 in V and -3.37 ± 0.06 
in K based on 58 stars. For the LMC they combine data for the 
OGLE Cepheid sample (see below) with data from Persson et 
al. (2004). They find a slope of -2.73± 0.03 in V, and -3.23± 
0.03 in K. The conclusion by F07 is that Galactic and LMC 
slopes are in good agreement. 

There are some issues however that may make a direct com- 
parison between the LMC and the Galaxy less obvious. The dis- 
tribution of the Galactic Cepheids that define the calibration is 



more or less uniform in log P (both in F07 and in the present 
study), while the majority of OGLE fundamental mode Cepheids 
has a log P between 0.4 and 0.7 and a tail to longer periods 
(Udalski et al. 1999). In addition it has been demonstrated that 
the PL-relation in the LMC is non-linear in V (e.g. Ngeow et al. 
2008 and references therein). Both effects make that the slope 
based on a least-square fit to the sample of calibrating Galactic 
Cepheids cannot be compared directly to the slope derived for 
the LMC sample. 

To investigate this further we made use of the data kindly 
made available by Dr Soszyiiski which consists of an OGLE 
sample of 701 fundamental mode Cepheids whereby the opti- 
cal light curves have been used to transform the single-epoch 
2MASS data into intensity-mean values following the prescrip- 
tion in Soszyiiski et al. (2005). This is a slightly improved ver- 
sion of the dataset that Dr. Soszyiiski made available to F07. 
The mean V magnitude and E{B — V) values are taken from 
the original OGLE data (Udalski et al. 1999). The 2MASS data 
were transformed to the SAAO system following Eq. 10 in Koen 
et al. (2007). Removing outliers with iterative 3<t clipping the 
slope in the LMC F-band PL-relation becomes -2.727 ±0.033 
based on 664 stars. 

In a Monte Carlo simulation a random value of log P was 
drawn between 0.45 (the lower limit of the Galactic Cepheids) 
and 1.6. From the OGLE sample the object closest in logP was 
selected and its V and K magnitude taken. In this way samples 
of 66 stars were generated, and a PL-relation fitted. This was 
repeated a large number of times. The mean slope and the dis- 
persion in the slope is —2.69 ± 0.06. 

In the i^-band, after removing outliers, the PL-relation 
based on 672 stars has a slope of —3.193 ± 0.027. Using the 
Monte Carlo method with samples of 66 stars the mean slope is 
-3.23 ±0.05. 

The conclusion is that, within the errors, the distribution of 
period and/or a non-linearity effect does not bias the derived 
slope in V and K. 

In V the observed slope in the LMC agrees within la with 
the one derived for the Galactic sample both for a constant p- 
factor and for the one with a mild period dependence. In K the 
agreement is less good, in the sense that the p-factor with a de- 
pendence on period agrees more closely. 

6. Summary and discussion 

The period-luminosity in the V and iiT-band and the de- 
pendence on metallicity are investigated for a sample of 68 
Galactic Cepheids with an individual metallicity determination 
from high-resolution spectroscopy. Distances are derived using 
the Baade-Wesselink technique using the most recent surface- 
brightness relation and estimates for the projection factor. 
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Fig. 3. Period-Luminosity relation in the V-band. The bottom 
panel shows the residual plotted versus metallicity. Lines indi- 
cate the fits to the data. 
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Fig. 4. Period-Luminosity relation in the ii'-band. The bottom 
panel shows the residual plotted versus metallicity. Lines indi- 
cate the fits to the data. 



The difference in slope between a constant p-factor (consis- 
tent with available interferometric data) and a mild dependence 
on period (following theoretical work by Nardetto et al. 2007) 
results in a difference in slope in the PL-relation which is of the 
same magnitude as the error in slope. 

The slope of the LMC Pi-relation in V and K is derived 
taking into account the difference in period distribution between 
the Galactic sample and the LMC sample. In V the slope for the 
LMC sample agrees with the galactic one for both assumptions 
regarding the dependence of p on period. In K the comparison 
of the slopes favours the mild period dependence of the p-factor 
on period. 

The metallicity dependence of the PL-relation is investi- 
gated and no significant dependence is found. A firm result is not 
possible as the range in metallicity spanned by the current sam- 
ple of galactic Cepheids is 0.3-0.4 dex, while previous work sug- 
gested a small dependence on metallicity only (typically —0.2 
mag/dex, see the introduction). 

It would be extremely valuable to analyse a sample of LMC 
and SMC Cepheids in a similar way. Unfortunately, the sam- 
ple of Cepheids with adequate optical, NIR and RV light curves 
(see e.g. Storm et al. 2004, 2005) for which a Baade-Wesselink 
distance has been derived (Gieren et al. 2005) and the sample 
of Magellanic Cloud Cepheids for which a metallicity has been 
derived from high-resolution spectroscopy (Mottini et al. 2008) 
have no overlap. High-resolution spectroscopy of MC Cepheids 
with Baade-Wesselink distances would therefore be extremely 
valuable. 
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